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This article concerns the Weyl series of spectral functions associated with the 
Dirichlet Laplacian in a d-dimensional domain with a smooth boundary. In the case 
of the heat kernel, Berry and Howls predicted the asymptotic form of the Weyl se- 
ries characterized by a set of parameters. Here, we concentrate on another spectral 
function, the (normalized) heat content. We show on several exactly solvable exam- 
ples that, for even d, the same asymptotic formula is valid with different values of 
the parameters. The considered domains are d-dimensional balls and two limiting 
cases of the elliptic domain with eccentricity e: A slightly deformed disk (e 0) 
and an extremely prolonged ellipse (e — > 1). These cases include 2D domains with 
circular symmetry and those with only one shortest periodic orbit for the classical 
billiard. We analyse also the heat content for the balls in odd dimensions d for 
which the asymptotic form of the Weyl series changes significantly. 
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1. Introduction 

We consider two types of spectral functions: The heat kernel and the (normalized) 
heat content, in the physical literature also known as the survival probability. They 
are associated with the spectrum of the Laplacian in a bounded domain. Let Q. 
be the domain in the rf-dimensional flat space, with a smooth boundary dO.. The 
spectrum of the Laplacian, say with the Dirichlet boundary condition (BC), is given 

by 

-A.^(r) = A.^(r) r e fi, 

<?!>(r) = r&dn. ^ ' 

The cigenvalucis < Ai < A2 < A3 • • • < Aj < • • • form a discrete set (Courant 
& Hilbert 1953). The associated eigenfunctions (pi,(f)2,(f'3, ■ ■ ■ ,4'j, ■ ■ ■ form an or- 
thonormalized basis of real functions, /j^ dr(^i(r)(/)j(r) = Sij. 

It is instructive to formulate the spectral problem of the Dirichlet Laplacian 
in the context of the diffusion (probability) theory. The conditional probability 
p{r,t\ro, 0) of finding a particle at a point r € at time t > 0, if it started from 
ro e at ^0 = 0, is governed by the diffusion (heat) equation 

^^^^^ = Ap(M|ro,0). (L2) 

This equation has to be supplemented by the initial condition p{r,t = 0|ro,0) = 
^(r — ro) and by the BC p(r,i|ro,0) = for r G dQ, which reflects absorption of 
the particle hitting the boundary. The conditional probability can be expressed in 
terms of the eigenvalues and eigenfunctions of the Dirichlet Laplacian as follows 

p{r,t\ro,0) = ^</>,(ro)0,(r)e-^^*. (1.3) 
3 

Various quantities can be constructed from the conditional probability. 

Historically, the most studied quantity was the heat kernel, defined as the trace 

K{t)= [ drdrop(r,f|ro,0)5(r-ro)= / dr p(r, t|r, 0) = V e"^^*. (1.4) 
Jn Jn J 

Another commonly studied quantity is the heat content, in the chemical physics 
called also the survival probability. The "local" heat content is defined as 

H{t;vo)= I drp(r,t|ro,0). (1.5) 
Jn 
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It represents the probability that the particle, localized at a point Fq G f2 at = 0, 
remains still diffusing in the domain fl at time t > 0, unabsorbed by the boundary. 

The BC is again the Dirichlct one: H{t, rg) = for Tq G dil. The normalized heat 
content is defined as the average of the local one over the whole domain, 



It represents the probability of finding the particle in 17 at time t > 0, if it was 
distributed uniformly with the density over the whole domain at to = 0. The 

normalization by l/|ri| allows us to study also infinite domains. In the mathematical 
literature, H{t) is interpreted as the amount of heat at time t inside the domain fl 
with the boundary dfl held at zero temperature, provided that at initial time the 
heat was distributed uniformly over fl. 

The heat kernel and the heat content have many similar properties. The most 
important results concern the small-f expansion of both quantities. The first three 
terms of these series are proportional to the volume \ fl\, surface \dQ\ and the inte- 
grated mean curvature, respectively. The heat kernel K{t) was intensively studied 
e. g. by Weyl (1911), Pleijel (1954), Kac (1966), Stewartson & Waechter (1971). 
The heat content was introduced by Birkhoff & Kotik (1954) who analysed its ID 
version. The first few coefficients of the small-f expansion of H{t) for Riemannian 
domains with both Dirichlet and Neumann BC were derived by van den Berg et 
al. (1993), van den Berg & Gilkey (1994) and DesJardins (1998). An important 
progress has been made by Savo (1998a,b) who derived a recurrence scheme for the 
coefficients of the small-t series. 

As concerns the heat kernel, it proves usehil to apply the Laplace transform 
with respect to t. For dimensions d > 1, the heat kernel has to be regularized 
by subtracting the leading terms of the f-expansion that would diverge under the 
Laplace integral. Thus we get the regularized resolventa K{s), see e.g. Stewartson 
& Waechter (1971), Berry & Howls (1994). For example, in the 2D case we have 

Kis)=iyte-'^[m-§^l. (1.7) 
The small-f expansion of K{t) corresponds to the large-s expansion of K{s): 

OQ (K) 

^(^)-E^- (1-8) 

n=l * 

This is the so-called Weyl series. 

It turns out that not only the small-n coefficients ci^^ contain the information 
about the domain characteristics. Based on a general Borel-transformed theory 
(Balian & Bloch 1972, Voros 1983), Berry & Howls (1994) developed a formalism 
for calculating the large-n coefficients of the Weyl series in 2D domains. They 
conjectured that 

+ „^oo, (1.9) 

where a and /3 are some parameters, I is the shortest periodic (stable) geodesies 
for a classical billiard in O. Note that the factorial nature of the coefficients makes 
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the Weyl series divergent. In the case of the disk domain of radius i?, for which the 
shortest periodic orbits with length I = 4i? form a continuous family, the authors 
predicted and tested numerically the parameter values a = l/{8y/2n) and /3 — 3/2. 
In the case of domains with only one (isolated) shortest periodic orbit, like the 
ellipse, the authors conjectured the increase of the parameter (3 by 1/2, i.e. /3 = 2. 

Later Howls & Trasler (1999) extended these results to higher dimensions. They 
derived exact a, /3 and the generalized interpretation of I in the case of d-dimensional 
balls of radius R. They had to distinguish between dimensions d according to their 
parity. For odd rf, the contribution of the shortest periodic orbit vanishes and the 
next shortest one, i.e. that with three bounces, becomes leading. Thus I becomes 
the perimeter of an equilateral triangle: 

The parameter /3 depends on d as follows 

r (5-d)/2 deven, 
^ " \ 7/2 -d dodd. ^ ' 

Explicit formulas for a{d) were also given by Howls & Trasler (1998,1999). Further, 
they proposed a generalization of the asymptotic formula (|1.9p which involves all 
periodic orbits Ij on the domain 51: 



j k=0 i 



Howls (2001) studied quantum balls threaded by a single magnetic flux at their 
center. While in 2D the parameters I and /3 are insensitive to the periodic orbits 
arising from the diffractive flux line, for a spherical domain these parameters are 
modified by diffractive orbits, in particular I = 2R and (3 — 3. 

No regularization is needed in the case of the Laplace transform of the heat 
content (van den Berg & Gilkey 1994). We can introduce the Laplace transforms 
for both the local heat content 

/"OO 

H{s;ro)= / dt e-'"' H{t,ro) (1.13) 

and the heat content 







His)= / dte-^''H{t)=J2 

Jo 7 = 1 



7? 



+ A,- 



(1.14) 



These Laplace transforms are related by 



^(s)-T^ / droH(,s;ro). (1.15) 



n 



The main goal of this paper is the analysis of the counterpart of the Weyl series 
(jl.Sp for the heat content. The Laplace transform of the heat content has the large-s 
expansion of the form 

oo 



n=2 ^ 
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We obtain the coefficients {c„} for few exactly solvable domains. For even dimen- 
sions d, the coefficients {c„} fulfill at asymptotically large n an analogy of the 
equation p.9p . 

Tin- 13 + 1) ,^ , 

Cn^a , n^oo, (1.17) 

where the power n — 2 in the denominator ensures that a remains dimensionless. 
The parameters a, (3 and I differ from those for the heat kernel. Our task is to 
determine these parameters for the studied domains and to point out their possible 
step- wise modifications under a symmetry change of the domain. A typical example 
of the symmetry change is a transition from a disk, possessing an infinite number 
of the shortest periodic orbits for the billiard, and an ellipse, with only one short- 
est periodic orbit. For odd dimensions d, the Gamma function in the asymptotic 
equation (jl.l7p has to be replaced by a bounded oscillating function (at least for 
the studied d-balls). 

In connection with Weyl series (jl.l6p for the heat content we mention the paper 
of van den Berg (2004) who showed that the relation between the coefficients c„ and 
the shortest periodic orbit does not hold in general. In particular, two domains with 
different shortest orbits can have the same c„ coefficients; the shorter of the two 
periodic orbits is determined by the difference of the two (exact) heat contents. This 
ambiguity of the Weyl series for the heat content inspires us to investigate the last 
for simple domains like d-dimensional balls and the ellipse. To our surprise, while 
a general theory seems to be more complicated for the heat content comparing to 
the heat kernel, the asymptotic Weyl series are explicitly available from the exact 
results and expansions for our simple domains. Similarities and differences between 
Weyl series for the heat content and the heat kernel are pointed out. 

The paper is organized as follows. In §2 we derive a differential equation for 
H{s; r). In §3 we write its exact solution and H(s) for d-balls; even dimensions are 
analysed in §3a and odd dimensions in §3b. §4 is devoted to the 2D ellipse domain 
of eccentricity e and the small- s expansion of H{s). In §5, for an arbitrary value of 
s, we solve exactly the leading terms for two special cases: A slightly deformed disk 
(e — )■ 0, §5a, §5b) and an extremely prolonged ellipse (e — ^ 1, §5c). We extract the 
large-n coefficients of the Weyl series from all exact solutions. §6 is the Conclusion. 



2. Differential equation 

The Laplace transform of the local heat content (|1.13p satisfies a differential equa- 
tion which is derived by the following procedure. Applying the conjugated Laplacian 
A+ (acting upon the coordinates of Tq) to the representation (|1.13p and using the 
conjugate of the diffusion equation (|1.2|) . we obtain 

A+#(s;ro)- / die"^'* / dr A+p(r, i|ro, 0) = / dr / dte-'"'dtp{r,t\ro,Q). 
Jo Jn Jn Jo 

(2.1) 

Integration by parts in t then implies the desired equation 

AH{s;t) - s'^H{s;r) = -1, (2.2) 

where we abandon the subscript and replace by A. This differential equation 
is supplemented by the Dirichlet BC H{s; r) = for r G d^l. A similar differential 
equation was derived by van den Berg & Gilkey (1994). 
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3. d-balls 



For d-balls, the Laplacian in (|2.2p can be expressed in terms of d-dimensional spher- 
ical coordinates. Since H^'''\s;r) does not depend on angle coordinates, we end up 
with the ordinary differential equation in r — |r|; see also van den Berg & Gilkey 
(1994). The solution reads 



d/2-l 



[sr) 



R'-'"^h,2-i{sR) 



(3.1) 



where I^{x) are modified Bessel functions. Integrating over the domain according 
to ([TTT5)) . we get 



1 



dld/2{u) 

uId/2-l{u) 



1 Id/2+l{u) 
Id/2^l{u) ' 



(3.2) 



where we introduced the scaled variable u = Rs. This result can also be found in 
the paper of van den Berg & Gilkey (1994), though in a slightly more complicated 
form. Its asymptotic analysis depends on the parity of d. 



(a) Asymptotic Weyl series for the halls in even d 



Our aim is to find the asymptotic large-n form of the coefficients c„ of the Weyl 
series (|1.16p resulting from the equation p.2[) . It is useful to rewrite the ratio of 
Bessel functions with the help of a logarithmic derivative: 



i?(^)(.) = l 



d{d-2) 

2w2s2 



d d 



From Gradshteyn & Ryzhik (2007) we use the large-u expansion 



Iu{u) 



r(j + i)r(i/-j + i) 



(3.3) 



(3.4) 



plus exponentially small terms; for the present v = c?/2 — 1 being integer, this is an 
infinite series. After straightforward calculation outlined in Appendix A, we get for 
asymptotically large n 



(-1) 



d/2-l_ 



Ad 



TT {2Ry 



r,„-3,+ (i^Mr(„-4, 



(3.5) 



Comparing with the asymptotic representation (|1.17p we see that the length 
I = 2R is now one half of the shortest periodic orbit l'^^^ — 4i? appearing in the 
case of the heat kernel. The symmetry parameter also changes significantly, /3 = 4 
regardless of (even) dimension. Contrary to the heat kernel, the whole dependence 
on dimension concentrates in the prefactor a^'^^ = (— l)'^/2~^4(i/7r. The next-to- 
leading term in p.Sp is analogous to the one in the heat- kernel series (|1.12D . 
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(b) Asymptotic Weyl series for the balls in odd d 

if^'')(s) in the formula p.2p with odd d involves Bessel functions /^(u) with 
half-integer f. The asymptotic series p.4p of such functions terminate, 

^" '^^"IP r(^ + J + ^) ^ ,,,, 

U(u) ^ , > -J—--- — : ^ — : TT, -u — > oo. (3.6) 

V2^^, {2uy r{j + i)r{u-j + ^)' 

For d = 1 we have 

HWis)^\(l--)+Oil/un- (3.7) 



Consequently, c„ = for n > 3. An analogous result holds for the Weyl series of 
the heat kernel. 

For d = 3, we also find a finite series 

H^'\s)^\(l-- + ^)+0{l/un, (3.8) 



so c„ = for n > 4. This differs from the heat kernel which has an infinite number 
of nonzero terms in Weyl series for all d > 2, see e. g. Howls & Trasler (1999). 
In the d = 5 case, the number of non-zero coefficients is infinite: 



J, <' - + '^"7 - '^"") , i, 1 - 1 + 12 _ 5f: 4], 

1 — U ^ S'^ \ U U'^ ^ — ' U*" / 

\ fc=3 / 

(3.9) 

so that c„ = — 5/i?" ^ for n > 5. We see that the asymptotic formula (|1.17l) no 
longer holds; There is no F-function (thus no /3) in the expression and the length 
I = R. 

For odd d > 5, the asymptotic formula for c„ acquires a new n-dependent term 
in the numerator. In contradiction to the function r(n — /? -f 1), this term exhibits 
bounded oscillations. The convergence of the series ^„ c„/s" will depend on R. 

For d = 7, these oscillations become periodic. Using the asymptotic expansion 
(|3.6p for /5/2(w) in the denominator of the representation p.2p and applying the 
partial fraction decomposition 

1 _ 1 / I ( ^ 

l-iv + iv"^ ^{vi - V2) \v - vi V-V2J ^ 3(wi - t;2) Vfi'^^ V2^^ 

(3.10) 

where v = l/u — l/{Rs) and vi — 3 exp(i7r/6), V2 — S exp(— i7r/6), we get 



R-" 



^ :3t-2 



fn'K\ r- /nTTXi 
cos J — v3sm J , n>5. (3.11 



The characteristic length is now I = R\vi\ = R/VS. The quasi-periodicity Cn+12 = 
3^R~^^Cn is due to the commensurability of the phase ±7r/6 of the complex conju- 
gate roots vi,2 with 27r. 

For d = 9, we have the polynomial 1 — 6v + 15v^ — 15v^ in the denominator. 
It has two complex conjugate roots vi, V2 and one real root V3, expressible in 
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Figure 1. Complex poles of H^'^\s) in the v — l/(Rs) plane 



terms of Cardano formulas. Numerically, |ui| = \v2\ ~ 0.39346201, with phases 
±0.2425136494068. . .TT, and V3 - 0.430628846. Since the complex phases are no 
more commensurate with 27r, the oscillations of c„ are not periodic. The reciprocal 
of the polynomial 1— 6w+15w^— 15w^ has three summands of the type w'^/w^^, m — 
1,2,3 as in the equation p.lOp . The coefficients |c„| with large n are dominated 
by the complex roots with the lowest absolute value |wi| — \v2\ < V3, thus we get 
I = R\vi\. A formal comparison with equation (|1.12p leads to the introduction of two 
lengths h = R\vi\ and I2 = Rv^, but we did not find any geometric interpretation 
of these lengths. 

For dimensions d > 9, we get more poles and more summands of type /v^^ . 
The leading term of |c„| is still given by the pair of conjugate poles with the lowest 
absolute value fi, W2 = v1 and the relation I = R\vi \ still holds. In figure 1, in the 
units of i? = 1, the roots of the denominator are plotted in the complex w-plane for 
larger dimensions d. The roots with negative imaginary part, placed symmetrically 
below the real axis, are not shown. The root vi with the lowest absolute value is 
always the leftmost point of the set. 

To conclude, the consideration of the large-s asymptotic for balls in odd d leads 
to H{s) which are the ratios of two polynomials and as such have a non-zero radius 
of convergence in s. On the other hand, the Berry and Howls conjecture (II. 9p as 
well as our conjecture (jl.l7p have zero radius of convergence in s and no longer 
hold. It is also questionable whether to expand the ratio of two polynomials into an 
infinite Weyl series, maybe the geometric information about the domain is hidden 
in the polynomial coefficients themselves. 

4. Ellipse, small-s expansion 

Now, in 2D, we shall pass from the disk to an ellipse which is an example of the 
domain possessing only one shortest periodic orbit. It will allow us to reveal a 
discontinuous change of the j3 parameter in (jl.l7D . 
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Let us consider the elliptic domain n centered at the origin, with major and mi- 
nor semiaxes a and b, respectively. In the Cartesian coordinates {x, y), its boundary 
is given by 

^ + § = 1, a>6. (4.1) 



The eccentricity is defined by e = ^Jl — {b/aY G (0, 1). The extreme cases e = 
and e — > 1 correspond to the disk and the infinitely prolonged (locally strip-like) 
ellipse, respectively. 

There is little hope to obtain explicitly the heat content for the ellipse domain 
with general e. Nevertheless, we are able to construct systematically several types 
of expansions for H{s;a,b). We derive the small-s expansion for arbitrary values 
of a, b in the present §4. Expansions with respect to e, around e = and 1, for 
arbitrary s will be analysed in the next §5. 

To find the small-s expansion of H{s) for the elliptic domain, we first formally 
expand the local H{s; r) in powers of s^, 

oo 

F(s;r)=^ijf (r)s2^ (4.2) 



Inserting this expansion into the differential equation (j2.2p implies an infinite se- 
quence of coupled equations obeyed by the unknown functions {ijj'*'' (r)}: 

AF(^'(r) = -1, (4.3a) 
AH!^'\r) = H^/\{r) for j > 1. (4.3b) 

Each of these functions must satisfy the Dirichlet BC 

H^/\r)=0 for re on, j==0,l,.... (4.4) 

It is convenient to work with complex coordinates z ~ x + iy and z — x ~ iy, in 
which the elliptic boundary (j4.ip becomes 

The Laplacian in complex coordinates has the form A — Adzdg- We solve succes- 
sively the set of equations (|4.3ap . (|4.3bp integrating their r.h.s. in z and z. Adding 
general solutions of the homogeneous equation A/ = 0, 

oo 

will permit us to fulfill the Dirichlet BC at the boundary (j4.5p . 
Starting with the equation (j4.3ap . we have 

Hi,'\z, z) = -\zz + + cW(,2 + ,2). (4.7) 



Article submitted to Royal Society 



10 



Igor Travenec and Ladislav Samaj 



The coefficients c^^ and Cq^'' follow from the condition H^^\z, z) = at the bound- 
ary (1331): 

The result (|4J)) with (|48| is substituted into the equation (j4.3bp for j = 1. After 
integrating and adding the homogeneous solution, we obtain 

+c^l'^+c['\z' + -z^) + c^;'\z^+-z'). (4.9) 

The coefficients cf \ cj^'' and c^^^ are again fixed to satisfy the Dirichlet BC for 

H'f \z,z) at the elliptic boundary. We proceed analogously in higher orders. The 
desired expansion of H{s) in powers of is obtained by averaging the relation (|4.2p 
over the ellipse surface: 

where = irab. The coefficients of the small-s expansion are obtained in the form 

^(3) _ {abf 
° " 4(a2 + 62) ^l.iij 

^1 - "l2(a2 + 62)2 (4-12) 

^ (ab)6[17(a^ + b^) + 98(ab)^] 

2 576(a2 + 62)3[(a4 + 54) + 6(a6)2]' ^ ' ' 
--(,) ^ (a6)S[93(a« + 6») + 1048a262(a4 + &4) ^ 3i90(nj,)4] 

3 ~ 8640(a2 + 62)4[(a4 + 64) + g(a5)2]2 ' ^ ■ > 

etc. 

For future purposes, we expand these coefficients around the limit e — )• 1, taken 
as a — 7- oo with b fixed. We include also the subleading term in 62/a2 = 1 — e2. 

^° " 4 4 U2 j ^ 



(,) _ 17 ,e 5566 ^^2^ ^^^^4 



2 576 576 la2 



3 2880 216 V a2 ' ' ' ^ ' 



5. Ellipse, eccentricity expansions 

We perform a change of variables x' = {b/a)x, y' — y; the Jacobian of this transfor- 
mation is J = a/b. The boundary then becomes x'^ + = b^, i- e. the transformed 
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domain fl' is the disk of radius b. The differential equation (|2.2[) modifies to 

H{s;r') - s^H{s;r') = -1. (5.1) 



dx' 



Within the probabihstic context explained in Introduction, instead of an isotropic 
diffusion in the "anisotropic" ellipse we get an anisotropic diffusion in the isotropic 
disk. This approach was inspired by the work Kalinay & Percus (2006). 
The differential equation (jS.ip can be formally expressed as follows 



(5.2) 



where A is a smallness parameter and A, B are the corresponding operators. There 
are two natural choices of the smallness parameter A. In §5a we shall set X — 
whereas in S5c we shall choose A = 1 — e^. In the case \ — e^, we have 



A + AS ) H{s\ r') - s^H{s; r') = -1, 



A = 



92 



92 



For A — 1 — we have 



92 



dx' 



A 



dy'- 



B 



dx' 



(5.3) 



(5.4) 



We look for the solution of equation (|5.2p perturbatively as an infinite series in the 
smallness parameter A: 



(5.5) 



n=0 



Inserting this expansion into (j5.2p and collecting terms of the same powers of A, we 
get a coupled set of differential equations 



AhI^\s-v')-s^H'^>{s-v') - -1, 
AHl^\s:v')~s'Hi^\s:v') = -BHi'\is;r') n = l,2 



2fr(£)/ 



(5.6) 
(5.7) 



All Hn \s;r') satisfy the Dirichlet boundary condition on the disk domain fl'. The 
quantity of interest H{s) is given by His) — X^i^o Hn\s)X"', where 



is-r'). 



(5.8) 



(a) Ellipse, small-e expansion 
We first choose A = as the smallness parameter. Equation (15.61) then becomes 
92 92 



dy'' dx'' 



H^°\s;r')-s'H^''>{s;r') = ^l, 



(5.9) 



where the upper index (0) refers to the e limit. The BC is i7^°^(s;r') = 
for |r'| = b. The solution in polar coordinates {r',(p') reads HQ^\s;r' , ip') = [1 — 
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/o(sr')//o(s6)]/s^ which is the special case of the equation p.l[) for d = 2 and 
R~h. Applying the relation (|5.8p and using the notation u = sb, we find the 2-ball 
version of the equation (j3.2p 



fv(o) _ 1 

^0 



1 



2Iiiu) 
uIo{u) 

The equation (|5.7p with n = 1 takes the form 



(5.10) 



i?(°)(.;r')-.^i?r(.;r') = -^i?r(.;r'). 



.2rr(0) 



(0)/ 



In polar coordinates, the r.h.s. becomes 



dx' 



/o(sr') /2(sr') 



+ e 



-2i(p' 



The general solution has the form 



f \s- r', = /io(s; r') + his; r') (e^^"^' + e"^*'^') . 



(5.11) 



(5.12) 



(5.13) 



hQ{s;r') is determined by the differential equation 

-^/io(s; r') + ^-^ho{s] r') - s^ft.o(s; r') 



2Io{sb)- 



(5.14) 



The homogeneous solutions are Io{sr') and KQ{sr'), their Wronskian is Vt^ = — l/r'. 
Using indefinite integrals like J z/o(2;)A'o(^)dz = 2;^[ifi(z)/i(z)+Xo(z)/o(z)]/2 that 
can be found in Bateman & Erdelyi (1953), the solution can be simplified to 



4s/o(s6) 



(5.15) 



We set = as we demand regular solution inside the domain, including the 
origin r' — 0. The constant c/ is fixed by the boundary condition /io(s; h) = 0. We 
thus get 



ho{s;r') = - 



r'hjsr') bh{sb)Io{sr') 
Aslaisb) 4s/2(s6) ' 



(5.16) 



Since the integral /q'^c^*''^ dip' — unless / — 0, only ho{s;r') contributes to 
and we finally obtain 



(0) 



fV(0) _ 1 
^1 -"2? 



2h{u) Il{u) 



uloiu) I^{u) 



(5.17) 
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Proceeding analogously in higher orders we find 



^ - 16.^.34)/|(.) + + 

m(24+ llu2)/3(y)/2(y) + (16 + 8^2 - 3w^)/^(u)/f(u) 
2u(2 + 5u2)/o(w)/^(w) - Bu^/f (u)] , (5.18) 
1 

192s2it3/4(w)/|(uj 
-2w(72 + 57^2 + lOu^)I^{u)lf{u) + 32(3 + Au^ + u'^)I^{u)lf{u) 
+m(8 + 52u^ + 15u*)/o2(m)/j'(m) -36m2(2 + u^)Ioiu)I^{u) + 2Au^I^{u)] 

(5.19) 

and so on. Recall that H{s) = J2T=o Hf\s)e^K 

We would like to emphasize that the obtained e-expansion of H{s) is valid for all 
values of s. This enables us to perform a consistency check of the above results by 
expanding them in small s and comparing with the previous small-s formulas (j4.1ip - 
(|4.14p . Expanding our Hf\s) in 1/s we can also test our results by comparison 
with the exact recurrence scheme of Savo (1998b) for the small-n coefficients c„. 
Our results pass also this consistency check. 

We can now analyse the large-s behavior of the set (I5.17p - (l5.19p . in analogy 
with Appendix A. Since analytic calculations are cumbersome, they were checked 
numerically as well. We found that the leading large-n term for the coefficient c„ 
coming from Hj^^ e'^'^ is proportional to e'^^[n'^ + 0{n'^~^)]. Collecting only these 
leading terms, we get 



c„(£) 



8r(n- 3) 
7r(26)«-2 



2\3 



4 64 768 



(5.20) 



Here, n ^ N where iV is a large number. So far, the series (|5.20l) is a formal 
expansion in ne^. As we are interested in the asymptotically large n for a fixed £, 
i.e. ne^ — >■ oo, we have to know all terms of the series (|5.20p . In the next subsection, 
we propose another method for finding H(s) based on plausible, but not rigorously 
justified, arguments. This method will predict all terms of the expansion (j5.20p . 
reproducing correctly the lowest ones. 



(&) Renormalized small-e expansion 

We return to the original (non-transformed) space r. The ellipse boundary R{ip) 
is expressed in polar coordinates as follows 

Riif) = , ■ f e (0, 27r). (5.21) 

y/1 — £^ COS^ (y9 

For small e, the ellipse is very close to the disk. Our intuitive approach is based 
on the assumption that in the differential equation (|2.2p for ll[s\ r) we are allowed 
to neglect the angular part of the Laplacian, i.e. A = 9^ + (l/r)9r- The resulting 
equation is equivalent to that of the disk, the dependence on the angle is included 
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only via the elliptic BC a.t r — R{(p): 



1 - 



loisRicp)) 



(5.22) 



In polar coordinates, the averaging over the ellipse domain is expressible as 



^/dr.... 



2 /-STr 



dip 



After the integration over r, we get 



H{s) 



2 /-Stt 



d(y5 



1 



dr r • • • 



2h{sR{ip)) 



(5.23) 



2s2(l -e2cos2 95) sfo^l -£2 cos2 lfIo{sR{ip)) 

' ... (^-^^^ 

To analyse the large-s behaviour of the ratio under integration, we repeat all steps 
(|A.2p - (|A.5p of Appendix A, the case = 0. The only difference compared to the 
2-ball consists in the replacement of the Bessel functions argument u by sR{ip): 



1 (1 - £2 C0S2 (py/'^ _ 1 

iP {sby ^ {sby 



^ ln(l - cos^ f) 



1 



-e 2 



(5.25) 

where we expanded the logarithm only up to the leading £2 term. The integration 
over the angle gives 



27r 



dip exp ( — cos^ 'p]= exp ( — ) Iq ( 



4 



(5.26) 



Similarly as in Appendix A, the l/s" term in the expansion of H{s) is identified 
with the substitution j = n — 4. Using that e~''Io{x + q) ^ laix) for large x, we 
finally arrive at 

T 



. ( X 8r(n-3) 2^ A 
^"^^^ ^(26)"-2 «^ ^0 



(5.27) 



Let us first assume that ne is finite. Then 



C„(£) 



8r(n-3) ^ (2fc- 1)!! / n£2" 



7r(26) 



n-2 .il^ ' 



fe=0 



(fc 



1^2 



(5.28) 



The first four terms of this expansion match perfectly those in the equation (|5.20p . 
In the limit of interest n£2 — > oo, according to p.4p it holds 



C„(£) 



8r(n-3) 2 16r(n-|) 
7r(26)«-2 y/2^£ ^ y/2^{2b)^^^e' 



(5.29) 



The last expression exploits the property r{n — 7/2) ^ r(n — 3)/y/n for n — !• oo. 
The comparison with the equation (11.17^ implies the parameter /3 = 9/2 for the 
ellipse with small e. This value is by 1/2 larger than /3 = 4 of a disk, in close analogy 
with the asymptotic Weyl series for the heat kernel. 
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For the heat kernel, Berry & Howls (1994) expected the parameters a and f3 to 
be of the order 1. Our result (I5.29P suggests something different for the heat content. 
We expect the divergence a{e) oc 1/e in the symmetry breaking limit e — !• 0+. This 
is an acceptable price for the step-wise change of (3 when restoring the circular 
symmetry. But this is a minor comment, the more important statement about the 
universality of /3 for domains with only one shortest periodic orbit still holds. 



(c) Ellipse with £ — > 1 

Now let us consider A = 1 — as the smallness parameter in the formalism 
developed at the beginning of §5. Equation (|5.6p takes the form 



dy 



^H^,'\s;r')-s'H^'>{s;r') = -l, 



(5.30) 



where the upper index (1) refers to the limit e — > 1. From the homogeneous solutions 
we choose only cosh(s?/'), as the odd function sinh(s?/') cannot fulfill the boundary 
condition at y' = ±\/W—x^. Thus we get 



1 - 



cosh(s?/') 



cosh 



,2 _ 2./2 



(5.31) 



Now we perform the averaging (|5.8p . considering for simplicity four times the first 
quadrant: 



da;' 



cosh(sy') 



cosh(s\/62 ~ x''^)_ 



(5.32) 



The integration over y' is simple. To integrate over x', we make the substitution 
x' = b sin Lp and resort to the full angle integration. 



2ir 



COS (p tanh(6s cos (p)d(p. 



(5.33) 



This is the exact Laplace transform of the heat content for the limiting case of an 
ellipse with finite width 2b and infinite length 2a — >■ oo, valid for any value of s. 

To check this expression within the small-s expansion, we apply the series rep- 
resentation (Gradshteyn & Ryzhik 2007) 



tanhz.^iii^i?,,z--. 



(5.34) 



where B2k are Bernoulli numbers. Inserting this series into (j5.33p . the integration 
results in 



(22'= - l)B2k 



k=2 



(fc!) 



1^2 



{bs 



,2k-~4 



(5.35) 



This series can be compared with the leading terms in the set (j4.15p and we find a 
perfect agreement in all available orders. It is worth mentioning that for complex s 
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the series (|5.35l) converges if |s| < 7r/26 VM, i- c. up to the first imaginary poles 
given by the lowest eigenvalue in (|1.14p . 

The next-to-leading term H[^\s,r') fulfills the equation (j5.7l) with n = 1: 



dy' 

With respect to (|5.3ip . the r.h.s. is equal to 



dx' 



(5.36) 



92 



■j^H^Q\s;r') = cosh(sy'; 



92 



1 



where we introduced A{x') for brevity. The solution of (j5.36p reads 



cosh{sy')A{x'), (5.37) 



£•(1) , A{x') . A{x') r— 
H{ (s;r ) = — — -ysmh(sy) ^ — -vb^ 



2s 



2s 



x' tanh 



(^sVb^ - x''^^ cosh(sy'). 



(5.38) 

To calculate h[^\s), we first integrate over y' in analogy with (j5.32l) , integrate by 
parts with respect to x' and substitute x' = bsimp, to get 



dip sin ip 



sinh2(&s cos(y9) 



sinh'^(6s cos (p) 



cosh^ {bs cos ip) bs cos ip cosh^ {bs cos ip) 

(5.39) 

To check this formula, we expand the integrated function in powers of s to get 



55 b^ 



-s^ 

6 576 



11 5^ f, 4487 
-s — 



216 



172800 



(5.40) 



The first four terms can be compared with the b'^ /a? terms in the set of four 
equations (|4.15p and we see the full agreement. 

Now we are ready to analyse the large-s expansion of the exact solutions (|5.33p 
and (|5.39p . The calculations are presented in Appendix B. Except for the obligatory 
l/s^ term, only the odd powers of 1/s appear in the Weyl series. The results for c„ 
with n odd are summarized by the formula 



i + \ + o{\') 



16 



2 r(n- 7/2) 



7r3 {2by 



n — )■ 00. 



(5.41) 



Comparing with the representation (jl.l7p . we see that I — 2b, i.e. one half of the 
shortest periodic orbit, as was expected from the small-e analysis. The symmetry 
parameter /3 = 9/2 is reproduced as well. The prefactor is non- universal, dependent 
on X = 1 ~ e'^. 



6. Conclusion 

This paper concerns the asymptotic form of the Weyl series for the heat content 
associated with the Dirichlet Laplacian in a smooth domain fl. Using the meth- 
ods developed by Balian & Bloch (1972) and Voros (1983), Berry & Howls (1994) 
mapped the quantum billiard model onto the resolventa of the heat kernel and con- 
jectured a "universal" geometric interpretation of the parameters / and /3 in high 
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orders of the Weyl series (11.91) for general domains. It is questionable whether an 
analogical approach to the heat content is possible. Some doubts come from the 
finding of van den Berg (2004) that two domains with different shortest periodic or- 
bits can have the same Weyl series for the heat content. One can imagine a scenario 
analogous to that for the heat kernel where unstable periodic orbits are excluded 
from the formalism. Maybe the accessibility conditions for orbits are even more 
restrictive for the heat content; they might be satisfied for even-dimensional balls 
and the ellipse, but no more for an annulus or twice-cut-disk from van den Berg 
(2004) examples. One can also imagine a general analysis of the asymptotic Weyl 
coefficients starting from Savo's recurrent scheme (Savo 1998b). 

Our strategy was to analyse the asymptotic Weyl series for the heat content, 
conjectured in the form (jl.l7p . from the exact results for simple domains. These re- 
sults were obtained by solving the differential equation (|2.2I) with Dirichlet BC. For 
balls in even dimensions d, the conjecture p.l7p applies when we identify I — 2R 
and /3 = 4 independently of d. For balls in odd dimensions, H{s) is the ratio of two 
polynomials and our conjecture no longer holds. It might be that the geometric in- 
formation about the domain is contained in the polynomial coefficients themselves. 
Another open problem is whether the symmetry of balls is more important than 
dimensionality, or vice versa, when adapting our results to non-ball domains. 

Further we studied the ellipse, which represents domains with single periodic 
orbit, in two limiting cases of eccentricity e — and e — >■ 1. In both cases, the 
parameter /? is shifted by 1/2 compared to the disk. This phenomenon is in close 
analogy with the heat kernel. 

We are grateful to P. Kalinay for valuable discussions. This work was supported by the 
Grants VEGA No. 2/0113/2009 and CE-SAS QUTE. 



Appendix A 

We are interested in the asymptotic terms of the Weyl series implied by the equation 
p.3|) . From (j3.4p . we rewrite the large-u asymptotic of I,y{u) with v = d/2 — 1 as 
follows 

'"^"^ 72^^' + ''^' ''-^^(2.).T(, + i)r(.-,-fi)- (^-'^ 

Since S is small in the large-u limit, we can expand 

\n{2TTu) + E - -S^ + • • • . (A.2) 

Further we will explore the identity r(l/2 — j)T{l/2 + j) = (— l)-'7r for integer j. 
In what follows, we shall argue that in the series (|A.2p E contributes to the leading 
and higher-order terms, E^ contributes to the sublcading and higher-order terms, 
etc., of the coefficients c„ in the asymptotic Weyl series. 

Let us first analyse just the leading term given solely by E^. Substituting the 
truncated expansion (jA.2|) into p.3p and using standard properties of F-functions, 
we get the leading contribution 

Hi^)(,) d ^ 1 (i^+j)...(| +,■)(! +,) T^l+j) 

^,2 Z.2.V+i(^-,)(^ -,)...(! -j) nVij) ■ ^ ■''> 
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Further we need the large-j expansion 

' , ^ ^ 1 - — + . . . J » 1. (A.4) 

ro-)r(j + i) 4j ^ ^ 

Taking the Hmit _;/—)■ oo of the long ratio in (jA.3|) we get ( — 1)''/^^^. Considering 
the leading (first) term of the expansion (jA.4[) . we obtain 

s..M.,_,,./.-._iLf_L_r(z±i). 

j=i 

Comparing this series with (|1.16l) . we set j = rt — 4 and find the leading term of 
the asymptotic formula p. 51) . 

The calculation of the subleading 1 /ri term in the equation p. 51) is more com- 
plicated and, for simplicity, we restrict ourselves to the d — 2 case. There are two 
contributions. The simpler one comes from the subleading term in (|A.4[) . i.e. we get 
the leading factor 8r(n - 3)/7r(2i?)"-2 times (-l/4n). The tricky part comes from 
the term in (|A.2I) which contributes to c„+2/[s^(2m)"] by the sum 



^ 27r2r(fc + l) r(n-fc + l) 

k—1 



^ /I 1 1 1 

2 — -iFo —, —, —n: n, n; —1 

'^^\^22 2 2 

27ri (1 + n) 

Here, we introduced the hypergeometric function 3F2, see e.g. Gradshteyn & Ryzhik 
(2007). Its large n asymptotic is 

3F2Q,i-n;i-n,i-n;-l)^2 + i- + o(-l), n ^ 00. (A.7) 

Inserting this into (jA.6|) . in equation (13. 3p we get exactly the same contribution as 
above, i.e. the leading factor times (— l/4n). Summing up two equal contributions 
we get finally 8r(n-3)/7r(2i?)"-2 times (-l/2n); noting that r(n-4) - T{n-3)/n 
for large n, this is already the subleading term of the equation p.5p for the special 
case d = 2. One can show that E"^, etc. contribute to higher-order terms. 
The derivation can be generalized to higher dimensions d. 

Appendix B 

We aim at analysing the asymptotic Weyl expansion of the equation (|5.33p . Let us 
return to the first quadrant integration and rewrite appropriately tanh. 

Further we apply the substitution z — 26s cos ip and subsequently the series 
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We need to calculate the integral Jg dz z^^^^ j (e^ + 1). For very large s, the upper 
integration limit can be extended to infinity, because the added term is exponen- 
tially small 

0< / ^—<\ ^= 2fc + l le-^^'' ^ B.3 

We can calculate the large-fc asymptotic of 

—^{2k-\)\ k^oo. (B.4) 

+ 1 

Further we use the obvious relation {2k — 1)!! ~ (2fc)!2^'^/fc!. Applying the above 
steps, we get for both s and k large 



Considering the asymptotic behaviour of 

r(2fc + i)r(2fc + 3) o.,7r,2fc 



22fer(A: + l)r(fc + 2) V TT 



l + O 



A; ^ cx) (B.6) 



and identifying n = 2fc + 5, we arrive at the first term in the equation (|5.4ip . Note 
that n is odd, the even powers of 1/s do not appear. 

Subleading term in X = 1 — 
Let us first rewrite the equation (j5.39p in the following way 

where w — bs cos (p. F{s) fulfills the relation 



^ [sFis)] ^ 3 r ^7l-(^)'^. (B. 
ds ^ ^ ^' Jo bs V \bsJ cosh% ^ 



We expand 



' «, n2 sinh^ (2fc-2)! 

^bs J cosh' yj A.22fe-ifc!(fc-l)! UJ JJ. ^li.Jj 



The upper limit of the integration in (jB.Sp can be again extended to infinity, as in 
the case of the equation (|B.3p . Then we calculate [si^(s)]' and after all we get 

^^^ (2fc~2)!(2fc-l)! 1 
^(^) - 22^-ifc!(fc-l)! (2b.r+i - (^-^"^ 
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We substitute this series into (IB.7|) and use the large-fc behavior of the ratio 



(2fc- 3)r(2fc- l)r(2/c) 

22'=-ir(fc + i)r(/c) 




fc-^oo. (B.ll) 



Thus we obtain 

We set n = 2fc + 3 and finaUy arrive at the second term in the equation (|5.4ip . Note 
that only odd powers of 1/s appear again. 
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